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Abstract 

We study the transformation of the statistical mechanics of N particles to the 
statistical mechanics of fields, that are the collective coordinates, describing the 
system. We give an explicit expression for the functional Fourier transform of the 
Jacobian of the transformation from particle to collective coordinate and derive the 
Fokker-Planck equation in terms of the collective coordinates. Simple 
approximations, leading to Debye-Huckel theory and to the hard sphere 
Percus-Yevick equation are discussed. 
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1 Introduction 



Collective coordinates have provided a powerful insight into statistical 
problems notably in plasma physics [1], liquid helium [2], electrolytes [3] and the 
general theory of liquids [4, 5]. The use of collective coordinates raises several 
interesting problems, which we consider in this paper. 

It is clear that when a classical system of identical particles is considered, all 
the physically relevant observables that depend only on the coordinates, are 
functional of the density 

P(f) = £5(f-R B ), (1) 

a 

where R a are the coordinates of the particles (more general observables that depend 

also on momenta are functional of the density and the current density). We consider a 
system of N particles enclosed in a cubic box of size L and periodic boundary 
conditions. The natural collective coordinates are the Fourier transforms, with q ^ , 
of the density 

P q =^Jp(?)e I<i? d 3 r, with q,=^-. (2) 

An obvious problem in using the p q 's , is that their number is infinite, while 

the number of degrees of freedom is 3N . This means that only 3N are independent 
and all the other collective coordinates can be written as functional of the 
independent 3N . It is clear, however, that it is best to treat the collective coordinates 
in a symmetric way, keeping the infinite number of collective coordinates. We will 
show in the following, how this can be done by introducing the Jacobian of the 
transformation from particle to collective coordinates. 

An additional difficulty comes with the introduction of hard cores. Experience 
has shown that approximation of a hard core by a soft core is inadequate for, as in all 
problems which are basically topological, the values of thermodynamic and transport 
coefficients is dominated by the softness and 'hardening' of a soft core always leads to 
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false results. A major advance in the handling of the hard core was made by Percus 
and Yevick [4] and we discuss this in our context in section (3). 

We will discuss first the Jacobian of the transformation from particle to 
collective coordinates and show that its functional Fourier transform has a rather 
simple form [6]. Next we show that the problem of an explicit form of the Jacobian of 
transformation can be avoided by deriving an exact Fokker-Planck equation in 
collective coordinates. 



2 The Jacobian 



An average of a quantity G that is a function of the density is given by 



fridR Y G{pJexp 

Y 


[- 2 ^ T Ew(R.-R,)] 


fridR y exp 

Y 







Introducing 8 functions, we can write the numerator N in (3) as well as the 
denominator in the form 



N 



= f2)p(?)ndR T ri5 



f 



-iq R„ 



G{pjexp 



2kT 



(4) 



where 2)p denotes functional integration over general real functions p(f ) with fixed 

p = Vn, is the Fourier transform of the inter-particle potential, and p is the 
average density. (Note that p appears because of our definition of Fourier 
components of the density. Note also that since p is fixed the inclusion of the self 
energy causes no problem.) Defining the Jacobian 



j{p}=fridR Y n5 

J Y q>0 



-iq-R a 



NT 



(5) 



we see that N can be written as a functional integral over p 
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N = }2)p(r)J{p}G{pJexp -_^-£w A p 



(6) 



with a similar equation for the denominator. 

Consider now the functional Fourier transform of the Jacobian 



${cp}=j"2) P (r)j{ P }exp 



q>0 



(7) 



using eq. (5), we find that 



${cp}=[2) P (F)ndR T ri5 

where 



exp 



-i£<P q P- 



q#0 



= [rM] N ,(8) 



r{cp}= }dRexp[-i(p(R)], (9) 

and 

Since N is extensively large, this can be written also in terms of a fugacity 



so we could use ^{(p,(i} instead of ^{(p,N} 



${<P#}=exp[uT{(p}] 



(12) 



Thus 



J{p, u] = 91 |2)(p exp [pr{cp}] exp [i£ cp qP _ q ] , 



(13) 
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where Vi is the proper normalization and the sum now is over all q. 

Pertubation theory made by expanding (7) or (13) in cp yields a simple approximation 

for J{p}. We will not dwell on this here because we will obtain those approximations 

from the dynamical equations to be considered in the following. 

3. Dynamics 

The dynamic problem is of obvious interest in its own right, but also offers a new way 
to approach the hard core. The original study of dynamics of Bohm & Pines produced 
plasma oscillation i.e. a property of p , but the problems we have been interested in 
are better expressed by noise methods of Langevin and Smoluchowski. Although 
current usage is to call these equations Fokker-Planck, they are not really F-P for this 
classic work on Brownian dynamic involves the velocities (the ancient paper of 
Chandrasekhar [7] is still a splendid introduction to this equation.) We want to use the 
p and cp directly and so base ourselves on the equation for P(Roc...,t) 



a i aw 

■+- 



3R kT dR 



P = 



(14) 



w 



whose equilibrium solution, of course, is the Gibbs distribution P«e kT 
Now introduce P{p} via 



p{ P ,t}= fndRn 

J a q*0 



iq-Rp 



P{R p ,t} 



(15) 



Therefore, multiply (14) by II 

q*0 



-iqR B 



and integrate all over R to obtain 



ap 
at 



D V(k-n-?-p Ap+yok 2 — Pk p 
Ntf ; ap7 k+J a Pj k dp. 



^Vn v j; ap7 k ~ jHj kT 



(16) 
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Similarly one can write this in Fourier variables cp 



+1 ^L( k ■ j)*P-k ^ — 9-jP + L Dk ^ — ^p 



at 



a w 



acp_ k+j a<p_j kT 



ip=o 



(17) 



and verify that the steady state solution is 



P °= exp 



I> ^ cp_ k + Nlogje- (R) dR 



4W 



(18) 



The main strength of the above equations is that it may be used to obtain directly 
hierarchies of equations for average quantities. Let F be a functional of the density 
that is not expected to diverge too strongly, like p q p_ q , P e P e P- e e etc. Suppose we 

ap 

are interested only in equilibrium properties. Then — = , we multiply eq. (16) by F 

at 

and integrate by parts to obtain 



2 1 aF \ D r/l J a 2 F 

£ Dk 2/ p £( k ) p ) + 

r \ap k / vnt? \ a p k pj 

pD r ,, ..wjaF \ n 



(19) 



It is clear that the p q 's are dependent (being an infinite set of coordinates in 

contrast to the 3N components of the R a 's). Therefore, an infinite number of 

relations must hold among the p q 's . This fact may lead to interesting sum rules. 

A similar equation can be obtained for the dynamical problem in which a 
quantity A is measured at equilibrium, the system is allowed to evolve freely and 
then a quantity B is measured at a later time t . We stick at present to the equilibrium 
problem. 
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If a 'random phase' approximation is made by keeping the sum over k and j only those 
terms where one of the p's is p = Vn (Note that derivatives are never with respect 
to p ) the equation collapses to 



— -VDk 2 — 

at r ap k 



dp 



a pw k p k 
+Pk + 



kT 



P = 



(20) 



with the equilibrium solution 



P eq oc exp 



z k 



1 + 



pWk 
kT 



PkP 



kH-k 



(21) 



that leads, for example, to the Debye-Huckel theory. 

The simple RPA form given above can be used to derive the Percus-Yevick 
[4] equation for hard spheres. To do this we have to think of the hard sphere system as 
an ideal gas with a constraint that forces the pair distribution function to vanish inside 
the hard core radius. The pair distribution function is given by 



g 2 ( ? i - A ) = =2"(p (?, )P (r 2 ) - pS (?, - A )) 

P 



(22) 



Imposing the constraint is done by adding to the Hamiltonian a term 
If/ 

— d 3 rdVX(r-r')p(r)p(r') where the Lagrange multiplier makes g 2 vanish for 
2 

|fj - r 2 |<R , the hard core radius. The Lagrange multiplier must vanish outside the 
hard core radius because the constraint is just within that radius. Denoting by X k the 
Fourier transform of X(f ) , we find that in eq. (20) W k is replaced by X k , so that 



PkP- 



1 



1 



1+ r^k 

kT 



P^k i+L 



(23) 
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where X k is the Fourier transform of a function that vanishes outside the hard core 
radius. Our equations are thus g 2 (r ) = for r < R and X (r) = for r > R [8], which 
are the Percus-Yevick equations for a system of hard spheres. 
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